Single-photon heat conduction in electrical circuits 
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We study photonic heat conduction between two resistors coupled weakly to a single superconducting mi- 
crowave cavity. At low enough temperature, the dominating part of the heat exchanged between the resistors is 
transmitted by single-photon excitations of the fundamental mode of the cavity. This manifestation of single- 
photon heat conduction should be experimentally observable with the current state of the art. Our scheme can 
possibly be utilized in remote interference-free temperature control of electric components and environment 
engineering for superconducting qubits coupled to cavities. 
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In nanoscale devices where the structure size is comparable 
to the characteristic wavelength of the heat and charge carri- 
ers, quantum mechanical effects can have a great impact on 
their transport properties. An important manifestation of this 
phenomenon is the quantized electrical conductance [1, 2], 
i.e., the maximum current that an ideal conducting channel 
can support is limited by the conductance quantum. Also the 
thermal conductance is quantized, a result which can be de- 
rived from very general principles concerning the limit of in- 
formation transfer [3, 4]. This general result implies that re- 
gardless of the nature of the carriers of thermal energy, there 
will be an upper bound on the amount of energy which can 
be transported by a single mode. This has been demonstrated 
experimentally for phonons [5] and electrons [6]. Recent ex- 
periments explicitly show that, as predicted in Ref. [7], the 
quantization of thermal conductance applies equally to pho- 
tons at low temperatures [8, 9]. Further proposals have been 
presented for harnessing this photonic channel in useful ap- 
plications [10] and a convenient circuit model to compute the 
transmitted semiclassical power has been introduced [11]. 

We extend these previous studies to the case of quantized 
heat transport mediated by single photons with a fixed fre- 
quency. This is achieved by using a superconducting cavity as 
the mediator for the heat. We derive the quantum mechanical 
interaction Hamiltonian between the cavity and the reservoirs, 
which allows us to describe the full quantum dynamics of the 
cavity modes under an engineered and tunable environment. 
The quantum model is found to agree extremely well with the 
semiclassical model for weak coupling and low temperatures. 
We note that the narrow bandwidth of the single-photon chan- 
nel allows for interference-free remote temperature control of 
electronic components working at different frequency bands. 

Microwave photons are extensively utiUzed in the var- 
ious implementations of superconducting quantum bits, 
qubits [12-15]. In these schemes microwave photons are em- 
ployed, for example, as carriers for transporting information, 
or in order to perform control operations on the qubits. Hence, 
microwave cavities form an integral component of these su- 
perconducting computer architectures and already impressive 
progress has been made on generating, guiding and manip- 
ulating individual microwave photons within electrical cir- 
cuits [16-19]. We demonstrate how a resistor embedded in 
such a cavity may be used as a controllable environment for 



the superconducting qubits. The control allowed by our tech- 
nique may therefore offer an effective method to gain further 
understanding and control [20, 21] of the decoherence pro- 
cesses of cavity qubits [22, 23]. 

The physical system we consider is shown schematically in 
Fig. 1(a). It consists of a superconducting coplanar waveguide 
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FIG. 1. (color online) (a) Schematic figure of two resistors, Ri and 
/?2> in a coplanar waveguide cavity, (b) The bare cavity can be mod- 
elled as a sequence of resistors and inductors with a capacitive con- 
nection to the ground. The values of r, (, and c are the resistance, 
inductance, and capacitance per unit length, respectively. This dis- 
crete model approaches the true continuous case as Ax — > 0. 

microwave resonator [24] with two resistors embedded close 
to the opposite ends of the cavity. We first consider the case 
of a bare cavity without resistors. Such a transmission line 
may be represented as a sequence of resistors, inductors, and 
capacitors as shown in Fig. 1(b). We neglect first the intrinsic 
dissipation of the cavity. The quantum mechanical voltage op- 
erator for this one-dimensional representation of a transmis- 
sion line can be found by quantizing the corresponding Euler- 
Lagrange equations in the continuum hmit [12, 25]. With the 
ends of the cavity defined at x = and x = L, and the boundary 
conditions imposed by charge neutrality, the resulting volt- 
age operator may be written as Vcav = 2^i^A'(0cos(fc7rx/L), 
where Ak{t) - yIfwJkKLc)[ak{t) + aj(f)], ^ is the length of the 
cavity and c is its capacitance per unit length. Here ak{f) and 
cifjii) are the bosonic annihilation and creation operators for 
the fcth mode, respectively. These operators satisfy the bosonic 

commutation relation [fl/t,a ] = 5ki and their time derivatives, 

J ^ 

in the Heisenberg picture employed here, are given by dtcik = 



-idJifCik and dio., - iciJka,, where a>k - kn/iLyfIc) is the fre- 
quency of the A:th mode in the cavity with { being the induc- 
tance per unit length. The current operator can be derived 
similarly and it is found to be /cav = Zi^i B kit) sin (kTix/L) -. 

ZZi h, with Mt) = / ^LhculKk'^n^Wkit) - alit)]. 

The Hamiltonian of the cavity-resistor system in the ab- 
sence of coupling is given by integrating the energy density 

in the cavity from to L as i7o = ^cav + ^^ - Jo ('^^cav + 
{I^^y)/2dx + Hr. We treat the coupling of the resistors to the 
cavity by adding their voltage fluctuations to the cavity volt- 
age operator in i/cav- Since the resistors are in practice orders 
of magnitude shorter than the photon wavelength they can be 
approximated as having a negligible spatial distribution and 
hence the potential drop across the ends of the resistor may be 
treated as a time-dependent discontinuity in the total voltage 
of the cavity at the position of the resistor If the resistor is 
placed at x', the Hamiltonian for a cavity coupled to an em- 
bedded resistor is 



H{t)^HR + - (ycav®Ires + Icav®5VLr^X 
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+ 2 I (Vcav®Ires + Icav®5VRr^X+- ( h^^dx. (1) 

Here 6Vufi) is the resistor-induced shift in the voltage on the 
left (right) side of the resistor and I is the identity operator. 
The tensor product, ®, has been included to emphasise that 
Vcav operates on the cavity and 6Vl,r on the resistor degrees of 
freedom. We neglect the fluctuations in the current due to the 
following reasons; The coupling between a resistor and the 
cavity modes decreases towards the ends of the cavity. Thus 
to achieve weak coupling the resistors are placed very close 
to the cavity edges. Since the current operator vanishes at the 
edges and the energy contribution from the current is propor- 
tional to square of the current operator, this contribution is 
much smaller than the one arising from the voltage. Hence, 
the energy fluctuations are almost entirely detemined by the 
fluctuations in the voltage. 

The Hamiltonian describing the interaction between the 
cavity and the resistor, Hint, is found by calculating the dif- 
ference in the energy of the cavity with and without the resis- 
tor. Neglecting second order terms in 6Vl,r yields the analytic 
result 



Ant = H{t) - Ho = - y -^dtlk ® <jt>res =: g ® 5%, 



(2) 



where (5Vres - 5Vl - SVr is the total voltage fluctuation across 
the resistor and Qix') = -2^1 (J^^dthix') is the integral of 
the charge density in the cavity from to .x'. Note that Hmx(t) 
is of the form Q ® 6V, where Q acts on the system and 6V 
on the Hilbert space of the environment. If the coupling 
strength is weak enough such that first order perturbation the- 
ory is valid and the autocorrelation time of the environment 
is sufficiently short as in a typical case, the transition rate 
between the cavity states \m) and |/) may be obtained from 



the Fermi golden rule as F,,,^; w \{l\Q\m)\^h~^Sv{-<j^mi) [26]. 
Here, w„,/ - (Ei - E„,)ITi corresponds to the energy change of 
the transition and S i/(w) is the spectral density of environ- 
mental fluctuations causing the transition. The voltage fluc- 
tuations can be described as a semiclassical voltage source 
with the Johnson-Nyquist spectral power density S sv^^^i^Jj) - 
2Rha)/[l -exp(pY)]. In the case of two resistors, with the 
first resistor at xi and the second at X2, the transition rates 
from the individual resistors can be added seperately since the 
coupling of the resistors to the cavity is linear and the volt- 
age fluctuations of the resistors are not intrinsically correlated. 
Thus the transition rates of the different photon states in the 
cavity assume the forms 
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where we have defined y*'^ = 2cjiR^^I^/(nZo) = 2R'^y(L{), and 

(i) 9 

an eff'ective resistance R ' - /?,sin (;7rx,/L). Here the char- 
acteristic impedance of the transmission line is Zq - y/TJc. 
Intrinsic dissipation mechanisms can also be included by giv- 
ing the cavity a small effective resistance per unit length r. 
This introduces a similar but far smaller term into the transi- 
tion rate corresponding to the finite Q value of the bare cavity. 
Note that y is inversely proportional to the total inductance of 
the cavity. Replacement of the center conductor with a series 
of superconducting quantum interference devices (SQUIDs) 
enables the tuning of £ independent of the other parameters. 
This would provide in situ control of the coupling strength and 
the transmission rates. 

We aim to find the quasiequilibrium electron temperature of 
the second resistor, T2, and the transmitted photonic power in 
the case for which the temperature of the first resistor, Ti, as 
well as the phonon bath temperature T, are fixed. Experimen- 
tally, it is possible to use superconductor-insulator-normal- 
metal (SIN) tunnel junctions to both control and measure the 
temperatures of the resistors [9, 27]. 

The master equation for the occupation probabilities 
Pm{t) - {p'"it)} of a cavity mode m can be expressed as 
LBmSH — Y'"p,„{t). The transition matrix has elements F"/; = 



dt 



ik 

Ff . for i + k and diagonal elements Ff, = -F'," , , , - F'," , , . 
The stationary state of this system corresponds to the zero 
eigenvalue of matrix F and yields the probability distribution 
in the quasiequiUbrium. 

To find the temperature of the second resistor we need to 
know only the net photon power transferred from the cavity 
to the second resistor, P\^^, and the power leaving the elec- 
tron cloud of the resistor to the phonon bath due to electron- 
phonon coupling, P\^ ^. Here, we assume that other heat con- 
duction mechanisms, such as the quasiparticle heat conduc- 



tion, are negligible. In equilibrium, there is no net heat trans- 
fer from the resistor and hence, P^j_ j^ = ^ cav- The electron- 

phonon coupling is given by ^^11,= £V(r^ - T^) where E is a 
material-dependent constant and V is the volume of the resis- 
tor [27]. Thus the quasiequilibrium temperature of the second 
resistor is given by 
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Since Pcav is the rate at which energy leaves the cavity due to 
absorption of photons by the second resistor, this is just equal 
to the energy of a single photon multiplied by the rate at which 
the number of photons in the cavity decreases due to Resistor 
2. For each mode, this quantity is the probability of having n 
photons in the mode multiplied by the energy difference in the 
transition and the net transition rate for that mode. To obtain 
the total power, we sum over all the modes as 
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In this paper F J' . represents the transition rate for the A:th 
mode due to the mth resistor. These are given by Eqs. (3) and 
(4) by keeping only the term corresponding to Resistor m. The 
first sum is over the modes and the second is over the photon 
number Since the transition rates depend on T2, Eqs. (5) and 
(6) must be solved self-consistently. The effective temperature 
of the cavity. Tea, can be defined in the spirit of the detailed 
balance as [26] 
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The parameters used are experimentally feasible: The ef- 
fective resistance and capacitance per unit length of the cav- 
ity are r = 2x 10"^ ilm"' and c - 130 x 10~^^Fm~^ respec- 
tively. The length of the cavity is L = 6.4 mm with cok - 
2nkx 10^" s~' . These values are consistent with the resonators 
fabricated in Ref. [24] and give a charcteristic impedance of 
Zq - 60.1 Q. Resistors 1 and 2 are identical with volume 
V - 2.25 X 10"^" m^ and are placed symmetrically in the cav- 
ity, ofi'set from the edge by 0.1 x L and 0.9 x L, respectively. 
We employ the parameters of Auo.75Pdo.25 with the electron- 
phonon coupling constant Z = 3 x 10^ Wm~^K~^. Similar re- 
sistors were produced in Ref. [9]. We compute the equilib- 
rium temperature for resistances of Ri = R2= R - 230 Q. and 
R^23Q. yielding y = 1.53 x lO** s"' and 1.53 x 10^ s"' re- 
spectively. 

The numerical results are shown in Fig. 2(a). At 40 mK 
bath temperature, Ti is changed from 40 mK to 400 mK. With 
R = 230Q. and at low Ti , the temperature of Resistor 2 closely 
follows that of Resistor 1 since the single-photon channel 
dominates over all other heat conduction mechanisms. The 
heat flow out from the hot electrons in Resistor 2 due to the 
electron-phonon coupling increases rapidly with the tempera- 
ture, and hence we observe a clear deviation of T2 from Ti 



above 250 mK. With R - 2.3 Q this heating eff'ect is still 
clearly observable although weaker due to the much weaker 
coupling. 

With a bath temperature of 250 mK, single-photon cooling 
of Resistor 2 is shown in Fig. 2(b). With R - 230 Q the tem- 
perature of Resistor 1 is lowered from 250 mK to 50 mK and 
T2 saturates to roughly 225 mK. At such bath temperatures 
the electron-phonon coupling makes a significant contribu- 
tion to the power transfer between the resistors and the bath. 
On the other hand, as the temperature of Resistor 1 decreases 
it excites exponentially fewer photons in the cavity and thus 
the photonic power exhange between the resistors saturates. 
As a consequence of the increased bath temperature, single- 
photon cooling is particularly susceptible to any reduction in 
the coupling, and hence is very low for R - 2.3 Q.. 

By altering either the resistance or the offset of the resis- 
tors from the edge of the cavity the coupling can be tuned. 
The stronger the coupling, the greater the change in the tem- 
perature of the Resistor 2, i.e., T2 follows Ti more closely at 
higher temperatures. However, our assumption of weak cou- 
pling breaks down beyond a certain limit, namely, the eff'ec- 
tive Q value of the cavity must be much greater than unity. 
For R = 230 Q, we have 2*=*^ ~ wi /(2y) = 20.5. 

We also analyzed the system semiclassically, employing the 
repeated lumped element model shown in Fig. 1(b) but with- 
out using quantum operators. The transmission line is divided 
into A^ points each with a capacitor of size cAx and an inductor 
of size £Ax. We take A^ = 100 which is sufficiently large for the 
discrete model to give a good approximation of the continuous 
transmission line. At two of these points k = {i,j}, the induc- 
tors are replaced with resistors and the corresponding voltage 
fluctuations are introduced. Applying Kirchoff 's current and 
voltage laws to the «th block in the circuit and eliminating 
the voltage gives, in Fourier space, the system of equations 
Z(w)/(w) = 6V{cij), where Z(aj) is a tridiagonal matrix which 
accounts for the impedance in the transmission line, 1(a)) and 
6V{cij) are vectors containing the current and voltage fluctua- 
tion at each point of the line. We hold Resistors 1 and 2 at 
constant temperatures and assume independent noise fluctua- 
tions. This allows us to calculate the change in the current due 
to the noise from each resistor seperately. The time-averaged 
noise power from Resistor / to Resistor j can be found directly 
from the current as 
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which can be calculated numerically for a given T, and Tj. It 
should be noted that Eq. (8) has the same form as the photonic 
power given in Ref. [7]. 

In the quantum case with hgT <sc Tiui, it is reasonable to 
approximate the cavity as a two-level system composing of 
zero or one photons in the fundamental mode. In this limit the 
photonic power transfer can be solved analytically if both T\ 
and T2 are known. The inset of Fig. 3 illustrates the process 
schematically. Since the probabilities must sum to unity, pQ + 
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FIG. 2. (color online) Temperature of Resistor 2 (dash-dotted line) 
and the effective temperature of the cavity (solid line) plotted against 
the temperature of Resistor 1 . The temperature of Resistor 1 (dashed 
line) is also shown for reference. The darker curves were obtained 
for R = 230 Q and the lighter curves for R = 2.3 Q.. In panel (a), 
single-photon heating is observed for a phonon bath temperature, 
T = 40 mK. In panel (b), single-photon cooling is observed for 
T = 250 mK. 
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Therefore the net 



quantum photon power transferred out of the second resistor 



Pq = r;_pohcoi - Y-^piTicji = ^ (rjr- - q r) , (9) 

which can be expressed in terms of the physical quantities us- 
ing Eqs. (3) and (4). 

The results of Eqs. (8) and (9) are plotted in Fig. 3 along- 
side the full numerical treatment of the quantum power The 
temperatures of both resistors are fixed here, with the tem- 
perature of Resistor 2 held 20 mK below the temperature of 
Resistor 1 which is increased from 40 mK to 150 mK. With 
R - 2.3 Q there is very close agreement between quantum and 
semiclassical power even as the power output changes over 
several orders of magnitude. At high resistances comparable 
to 230 Q, the quantum and semiclassical models have a some- 
what different temperature dependence. 
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FIG. 3. (color online) Quantum (broken lines) and semiclassical 
power (solid line) from Resistor 1 to Resistor 2 for R = 2.3 il (bot- 
tom lines) and R = 230 il (top lines). The temperature of the first 
resistor spans the range from 40 mK to 150 mK while the second 
resistor is held at at T2 = T\- 20 mK. The full numerical simula- 
tion (dashed line) which includes all the modes is in good agreement 
with the two-level approximation (dash-dotted line) in this tempera- 
ture regime. Inset: illustration of the two-level model for the system 
of Fig. 1(a). 



In summary, we have presented analytical and numerical 
results which show that heat conduction by single photons 
has a distinct signature which is readily observable in the heat 
transport between two resistors coupled by a superconduct- 
ing cavity. Photonic heat conduction enables the possibility 
to perform remote heating or cooling of low-temperature cir- 
cuit elements at an extremely narrow bandwidth, thus having 
a minimal perturbation to the usual operation of the compo- 
nent. Furthermore, we showed that although weak, the cou- 
pUng of resistors to the superconducting cavity can dominate 
the dissipation mechanisms, and hence completely determine 
the temperature of the cavity, i.e., these resistors work as en- 
gineered and controllable environments for the cavities. Thus 
our work shows that SIN thermometry and temperature con- 
trol can be extended to superconducting cavities. We observed 
a disagreement between the quantum and semiclassical mod- 
els for relatively strong coupling. This calls for further theo- 
retical and experimental studies on the true quantum nature of 
single-photon heat conduction. We thank Jukka Pekola and 
Teemu Ojanen for helpful discussions, and Academy of Fin- 
land and Emil Aaltonen Foundation for financial support. 
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